GeneWaltz--A new method for reducing the false positives of gene finding

<p>Abstract</p> <p>Background</p> <p>Identifying protein-coding regions in genomic sequences is an essential step in genome analysis. It is well known that the proportion of false positives among genes predicted by current methods is high, especially when the exons are short. These false positives are problematic because they waste time and resources of experimental studies.</p> <p>Methods</p> <p>We developed GeneWaltz, a new filtering method that reduces the risk of false positives in gene finding. GeneWaltz utilizes a codon-to-codon substitution matrix that was constructed by comparing protein-coding regions from orthologous gene pairs between mouse and human genomes. Using this matrix, a scoring scheme was developed; it assigned higher scores to coding regions and lower scores to non-coding regions. The regions with high scores were considered candidate coding regions. One-dimensional Karlin-Altschul statistics was used to test the significance of the coding regions identified by GeneWaltz.</p> <p>Results</p> <p>The proportion of false positives among genes predicted by GENSCAN and Twinscan were high, especially when the exons were short. GeneWaltz significantly reduced the ratio of false positives to all positives predicted by GENSCAN and Twinscan, especially when the exons were short.</p> <p>Conclusions</p> <p>GeneWaltz will be helpful in experimental genomic studies. GeneWaltz binaries and the matrix are available online at <url>http://en.sourceforge.jp/projects/genewaltz/</url>.</p

Optimality of mutation and selection in germinal centers

The population dynamics theory of B cells in a typical germinal center could play an important role in revealing how affinity maturation is achieved. However, the existing models encountered some conflicts with experiments. To resolve these conflicts, we present a coarse-grained model to calculate the B cell population development in affinity maturation, which allows a comprehensive analysis of its parameter space to look for optimal values of mutation rate, selection strength, and initial antibody-antigen binding level that maximize the affinity improvement. With these optimized parameters, the model is compatible with the experimental observations such as the ~100-fold affinity improvements, the number of mutations, the hypermutation rate, and the "all or none" phenomenon. Moreover, we study the reasons behind the optimal parameters. The optimal mutation rate, in agreement with the hypermutation rate in vivo, results from a tradeoff between accumulating enough beneficial mutations and avoiding too many deleterious or lethal mutations. The optimal selection strength evolves as a balance between the need for affinity improvement and the requirement to pass the population bottleneck. These findings point to the conclusion that germinal centers have been optimized by evolution to generate strong affinity antibodies effectively and rapidly. In addition, we study the enhancement of affinity improvement due to B cell migration between germinal centers. These results could enhance our understandings to the functions of germinal centers.Comment: 5 figures in main text, and 4 figures in Supplementary Informatio

Developing and applying heterogeneous phylogenetic models with XRate

Modeling sequence evolution on phylogenetic trees is a useful technique in computational biology. Especially powerful are models which take account of the heterogeneous nature of sequence evolution according to the "grammar" of the encoded gene features. However, beyond a modest level of model complexity, manual coding of models becomes prohibitively labor-intensive. We demonstrate, via a set of case studies, the new built-in model-prototyping capabilities of XRate (macros and Scheme extensions). These features allow rapid implementation of phylogenetic models which would have previously been far more labor-intensive. XRate's new capabilities for lineage-specific models, ancestral sequence reconstruction, and improved annotation output are also discussed. XRate's flexible model-specification capabilities and computational efficiency make it well-suited to developing and prototyping phylogenetic grammar models. XRate is available as part of the DART software package: http://biowiki.org/DART .Comment: 34 pages, 3 figures, glossary of XRate model terminolog

Chiral tunneling and the Klein paradox in graphene

The so-called Klein paradox - unimpeded penetration of relativistic particles through high and wide potential barriers - is one of the most exotic and counterintuitive consequences of quantum electrodynamics (QED). The phenomenon is discussed in many contexts in particle, nuclear and astro- physics but direct tests of the Klein paradox using elementary particles have so far proved impossible. Here we show that the effect can be tested in a conceptually simple condensed-matter experiment by using electrostatic barriers in single- and bi-layer graphene. Due to the chiral nature of their quasiparticles, quantum tunneling in these materials becomes highly anisotropic, qualitatively different from the case of normal, nonrelativistic electrons. Massless Dirac fermions in graphene allow a close realization of Klein's gedanken experiment whereas massive chiral fermions in bilayer graphene offer an interesting complementary system that elucidates the basic physics involved.Comment: 15 pages, 4 figure

An evaluation of enteral nutrition practices and nutritional provision in children during the entire length of stay in critical care

&lt;b&gt;Background&lt;/b&gt; Provision of optimal nutrition in children in critical care is often challenging. This study evaluated exclusive enteral nutrition (EN) provision practices and explored predictors of energy intake and delay of EN advancement in critically ill children.&lt;p&gt;&lt;/p&gt; &lt;b&gt;Methods&lt;/b&gt; Data on intake and EN practices were collected on a daily basis and compared against predefined targets and dietary reference values in a paediatric intensive care unit. Factors associated with intake and advancement of EN were explored.&lt;p&gt;&lt;/p&gt; &lt;b&gt;Results&lt;/b&gt; Data were collected from 130 patients and 887 nutritional support days (NSDs). Delay to initiate EN was longer in patients from both the General Surgical and congenital heart defect (CHD) Surgical groups [Median (IQR); CHD Surgical group: 20.3 (16.4) vs General Surgical group: 11.4 (53.5) vs Medical group: 6.5 (10.9) hours; p &lt;= 0.001]. Daily fasting time per patient was significantly longer in patients from the General Surgical and CHD Surgical groups than those from the Medical group [% of 24 h, Median (IQR); CHD Surgical group: 24.0 (29.2) vs General Surgical group: 41.7 (66.7) vs Medical group: 9.4 (21.9); p &lt;= 0.001]. A lower proportion of fluids was delivered as EN per patient (45% vs 73%) or per NSD (56% vs 73%) in those from the CHD Surgical group compared with those with medical conditions. Protein and energy requirements were achieved in 38% and 33% of the NSDs. In a substantial proportion of NSDs, minimum micronutrient recommendations were not met particularly in those patients from the CHD Surgical group. A higher delivery of fluid requirements (p &lt; 0.05) and a greater proportion of these delivered as EN (p &lt; 0.001) were associated with median energy intake during stay and delay of EN advancement. Fasting (31%), fluid restriction (39%) for clinical reasons, procedures requiring feed cessation and establishing EN (22%) were the most common reasons why target energy requirements were not met.&lt;p&gt;&lt;/p&gt; &lt;b&gt;Conclusions&lt;/b&gt; Provision of optimal EN support remains challenging and varies during hospitalisation and among patients. Delivery of EN should be prioritized over other "non-nutritional" fluids whenever this is possible.&lt;p&gt;&lt;/p&gt

smyRNA: A Novel Ab Initio ncRNA Gene Finder

Background: Non-coding RNAs (ncRNAs) have important functional roles in the cell: for example, they regulate gene expression by means of establishing stable joint structures with target mRNAs via complementary sequence motifs. Sequence motifs are also important determinants of the structure of ncRNAs. Although ncRNAs are abundant, discovering novel ncRNAs on genome sequences has proven to be a hard task; in particular past attempts for ab initio ncRNA search mostly failed with the exception of tools that can identify micro RNAs. Methodology/Principal Findings: We present a very general ab initio ncRNA gene finder that exploits differential distributions of sequence motifs between ncRNAs and background genome sequences. Conclusions/Significance: Our method, once trained on a set of ncRNAs from a given species, can be applied to a genome sequences of other organisms to find not only ncRNAs homologous to those in the training set but also others that potentially belong to novel (and perhaps unknown) ncRNA families. Availability

Topological Anisotropy of Stone-Wales Waves in Graphenic Fragments

Stone-Wales operators interchange four adjacent hexagons with two pentagon-heptagon 5|7 pairs that, graphically, may be iteratively propagated in the graphene layer, originating a new interesting structural defect called here Stone-Wales wave. By minimization, the Wiener index topological invariant evidences a marked anisotropy of the Stone-Wales defects that, topologically, are in fact preferably generated and propagated along the diagonal of the graphenic fragments, including carbon nanotubes and graphene nanoribbons. This peculiar edge-effect is shown in this paper having a predominant topological origin, leaving to future experimental investigations the task of verifying the occurrence in nature of wave-like defects similar to the ones proposed here. Graph-theoretical tools used in this paper for the generation and the propagation of the Stone-Wales defects waves are applicable to investigate isomeric modifications of chemical structures with various dimensionality like fullerenes, nanotubes, graphenic layers, schwarzites, zeolites

Current perspectives on invasive nontyphoidal Salmonella disease.

PURPOSE OF REVIEW: We searched PubMed for scientific literature published in the past 2 years for relevant information regarding the burden of invasive nontyphoidal Salmonella disease and host factors associated with nontyphoidal Salmonella infection and discuss current knowledge on vaccine development. The following search terms were used: Salmonella, non typhoidal/nontyphoidal, NTS, disease, bloodstream infection, invasive, sepsis/septicaemia/septicemia, bacteraemia/bacteremia, gastroenteritis, incidence, prevalence, morbidity, mortality, case fatality, host/risk factor, vaccination, and prevention/control. RECENT FINDINGS: Estimates of the global invasive nontyphoidal Salmonella disease burden have been recently updated; additional data from Africa, Asia, and Latin America are now available. New data bridge various knowledge gaps, particularly with respect to host risk factors and the geographical distribution of iNTS serovars. It has also been observed that Salmonella Typhimurium sequence type 313 is emergent in several African countries. Available data suggest that genetic variation in the sequence type 313 strain has led to increased pathogenicity and human host adaptation. A bivalent efficacious vaccine, targeting Salmonella serovars Typhimurium and Enteritidis, would significantly lower the disease burden in high-risk populations. SUMMARY: The mobilization of surveillance networks, especially in Asia and Latin America, may provide missing data regarding the invasive nontyphoidal Salmonella disease burden and their corresponding antimicrobial susceptibility profiles. Efforts and resources should be directed toward invasive nontyphoidal Salmonella disease vaccine development

Limited effect of patient and disease characteristics on compliance with hospital antimicrobial guidelines

Objective: Physicians frequently deviate from guidelines that promote prudent use of antimicrobials. We explored to what extent patient and disease characteristics were associated with compliance with guideline recommendations for three common infections. Methods: In a 1-year prospective observational study, 1,125 antimicrobial prescriptions were analysed for compliance with university hospital guidelines. Results: Compliance varied significantly between and within the groups of infections studied. Compliance was much higher for lower respiratory tract infections (LRTIs; 79%) than for sepsis (53%) and urinary tract infections (UTIs; 40%). Only predisposing illnesses and active malignancies were associated with more compliant prescribing, whereas alcohol/ intravenous drug abuse and serum creatinine levels > 130 mu mol/l were associated with less compliant prescribing. Availability of culture results had no impact on compliance with guidelines for sepsis but was associated with more compliance in UTIs and less in LRTIs. Narrowing initial broad-spectrum antimicrobial therapy to cultured pathogens was seldom practised. Most noncompliant prescribing concerned a too broad spectrum of activity when compared with guideline-recommended therapy. Conclusion: Patient characteristics had only a limited impact on compliant prescribing for a variety of reasons. Physicians seemed to practise defensive prescribing behaviour, favouring treatment success in current patients over loss of effectiveness due to resistance in future patients

Correlation between nucleotide composition and folding energy of coding sequences with special attention to wobble bases

Background: The secondary structure and complexity of mRNA influences its accessibility to regulatory molecules (proteins, micro-RNAs), its stability and its level of expression. The mobile elements of the RNA sequence, the wobble bases, are expected to regulate the formation of structures encompassing coding sequences. Results: The sequence/folding energy (FE) relationship was studied by statistical, bioinformatic methods in 90 CDS containing 26,370 codons. I found that the FE (dG) associated with coding sequences is significant and negative (407 kcal/1000 bases, mean +/- S.E.M.) indicating that these sequences are able to form structures. However, the FE has only a small free component, less than 10% of the total. The contribution of the 1st and 3rd codon bases to the FE is larger than the contribution of the 2nd (central) bases. It is possible to achieve a ~ 4-fold change in FE by altering the wobble bases in synonymous codons. The sequence/FE relationship can be described with a simple algorithm, and the total FE can be predicted solely from the sequence composition of the nucleic acid. The contributions of different synonymous codons to the FE are additive and one codon cannot replace another. The accumulated contributions of synonymous codons of an amino acid to the total folding energy of an mRNA is strongly correlated to the relative amount of that amino acid in the translated protein. Conclusion: Synonymous codons are not interchangable with regard to their role in determining the mRNA FE and the relative amounts of amino acids in the translated protein, even if they are indistinguishable in respect of amino acid coding.Comment: 14 pages including 6 figures and 1 tabl